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We present an approach based on the dynamical mean field theory which is able to give the 
excitation spectrum of a triply degenerate Hubbard model with a Hund's exchange invariant under 
spin rotation. The lattice problem can be mapped onto a local Anderson model containing 64 
local eigenstates. This local problem is solved by a generalized non-crossing approximation. The 
influence of Hund's coupling J is examined in detail for metallic states close to the metal insulator 
transition. The band-fllling is shown to play a crucial role concerning the effect of J on the low 
energy dynamics. 
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I. INTRODUCTION 

The Mott metal insulator transition (MIT) has been 
the subject of many experimental and theoretical stud- 
ies over the past decades.^'^ Among correlated materials 
where the Coulomb repulsion between electrons is strong 
enough to overcome their kinetic energy, we find many 
3d transition-metal oxides (see ref. 0) and alkali-doped 
fuUeridesi^ 

The infiuence of orbital degeneracy on the properties of 
strongly correlated electronic systems has been the sub- 
ject of numerous recent studies. It has been shown un- 
ambiguously that the orbital degrees of freedom are of 
primary importance for understanding the nature of the 
Mott MITpSiiaj^Jiii^ More recently, the interesting 
and controversial possibility of an orbital selective MIT in 
a multiorbital system has been investigatediiSiSiiSSiSiSi 
Hund's coupling plays an important role in these sys- 
tems but reliable dynamical quantities are not available 
for a spin and orbital rotational invariant Hund's cou- 
pling. Very recently a new approach using an extension 
of Wilson's numerical renormalization group (NRG) has 
been proposed:^ The exchange part of the Hund's cou- 
pling influences dramatically the physics of the two-band 
Hubbard model and changes the nature of the Mott tran- 
sition. 

In the present paper, we introduce into the triply de- 
generate orbital model a non zero Hund's coupling pa- 
rameter J and we take into account its full expression 
including the so-called spin-flip term. The only neglected 
term is the orbital flip contribution or pair-hopping term 
between different orbitals.^^ 

Even the simplest lattice model for correlated elec- 
trons, the one-band Hubbard model presents strong diffi- 
culties and approximations are necessary to solve it. The 
dynamical mean held theory (DMFT)^ provides a so- 
lution in the limit of large spatial dimensionjS Within 
the DMFT framework, a correct description of the Mott- 
Hubbard metal insulator transition is possible for the 
one-band Hubbard model.— 



In a previous work we have proposed an approach 
based on the DMFT which was able to deal with the dou- 
bly degenerate Hubbard model^^ with a vanishing Hund's 
coupling. Metal insulator transitions for integer total oc- 
cupations per site n = 1, n = 2 and n = 3 were observed, 
with different values for the critical parameters on the 
transition lines. Besides, the transfer of spectral weight 
between different bands was interpreted in terms of mi- 
croscopic dynamical processes. Here, in order to give a 
more realistic description of the spectral function of a cor- 
related metal with several orbital degrees of freedom close 
to the MIT, we propose a generalization of our previous 
approach. Also, we develop a systematic description of 
the effects of degeneracy, of the Hubbard interaction and 
the Hund's coupling on the triply degenerate model. 

The triply degenerate model is relevant for transition- 
metal oxides or other transition-metal compounds where 
the cubic crystal fleld splits the 3d level into the eg and 
the t2g manifolds whenever it is the t2g manifold which is 
partially occupied. An example is VOj, (see ref. |3^ and 
references therein). It shows a transition from metal- 
lic to semiconducting behavior for x close to 1 and the 

V2+ 

ion provides three electrons to the t2g orbitals. In 
that special case, however, not only correlation but also 
disorder effects seem to be very important, even for the 
stoichiometric samples. Other examples are the heavy 
fermion compound LiV204 which remains metallic up to 
very low temperature^S*^ and LiTi204 which is super- 
conducting below T = 13.7 Ki^ Our approach is also 
relevant for MgTi204 (see ref. Isol and references therein) 
which is metallic above T = 260 K. 

In our study we investigate the t2g-derived spectral 
function. We find three different effects in the spec- 
tral function: (i) band effects (which are usually already 
well described by the local spin density approximation 
(LSDA)), (ii) the Hubbard satellites (at an energy scale 
of U, already well treated by DMFT of the one-band 
model), and (iii) multi-peak effects connected with the 
Hund's coupling (energy scale J). It is the last aspect 
which is new in our approach in comparison to previous 
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studies. Keeping in mind that these multi-peak effects 
should be important to interpret realistic spectra, we an- 
alyze nevertheless a model situation in the present work 
with a semi-elliptic density of states for the uncorrelated 
problem. 

In the following section the model Hamiltonian is pre- 
sented, it takes into account explicitly the effect of Hund's 
coupling. This section also presents the application of the 
DMFT on this Hamiltonian and a method of resolution 
based on the non-crossing approximation (NCA). In the 
third section, the main results concerning spectral prop- 
erties are discussed. 



II. HAMILTONIAN AND SOLVING METHOD 



The triply degenerate Hubbard Hamiltonian, including 
a non zero Hund's coupling parameter J can be written 
(see for instance ref. HOT) as follows : 



where 
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The sum < i,j > runs over nearest neighbor sites of a 
Bethe lattice and a,b = 1,2,3 are the band indices, c^^ 
(respectively Ciaa) denotes the creation (respectively an- 
nihilation) operator of an electron at the lattice site i 
with spin a and orbital index a and riiao- is the occupa- 
tion number operator for spin a and orbital a. The hop- 
ping between different orbitals will not be investigated 
here since tfj* = —tSab- The ref. 12 presents a general 
Hamiltonian for the multi orbital Hubbard model. The 
only term neglected in our approach is the pair-hopping 
term between different orbitals. The prefactors in Eq.(^ 
ensure the Hamiltonian is invariant under any linear com- 
bination of the t2g orbitals. 

Within the framework of the DMFT, by integrating 
out all fermionic degrees of freedom except those for a 
central site i = o, the lattice model ^ can be mapped 
onto an effective impurity model. The corresponding 
Hamiltonian iJeff contains a local part ifioc and a part 
corresponding to the coupling to the effective medium 
that has to be determined self-consistently. We then 
write : 



H, 



off 



loc 
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and the coupling with the effective medium is 



kaa 

+ y b r . 



WS represents the hybridization between the site i = o 
and the effective medium corresponding to orbital a. ei 
is the band energy of the same effective medium. 6^- 

aka 

(respectively b^^^) is the creation (respectively annihi- 
lation) operator of an electron in the effective medium 
a. As in the doubly degenerate situation,^- the effective 
medium is a multi-component one that can be character- 
ized by the effective dynamical hybridizations : 



^"M = E 



A; 



k 



By writing the equation of motion of the Hamiltonian 
Hcs we obtain : 

Besides, on a Bethe lattice, the Green's function presents 
the following property 

Therefore, the self-consistent equations of the DMFT can 
be simply written like 

JdL0)=eGaa{^) ■ 

The local problem that has to be solved self-consistently 
is then formally similar to the doubly degenerate case 
previously investigated. However we have to deal here 
with a large number Ng of local impurity states \Si) — 
jai, Q!2, Oiz) where aa represents the electronic occupation 
(0, t, i or tl) of orbital a at the impurity site o. Indeed, 
the third orbital a = 3 increases Ng from 16 to 64 = 4'^. 
Another difficulty comes from the spin-flip term : due to 
this term, the basis of local states A — {\Si)} does no 
longer diagonalize ifioc- 

In the following, the 64 states that diagonalize i/ioc 
wiU be noted \Si). They form the basis A = {\Si)}. Tak- 
ing into account the degeneracy of local states \Si), it is 
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possible to decrease the basis dimension to 17, the corre- 
sponding states will be noted |to). Indeed the 64 states 
can be classified in 13 generic families (third column in 
TableU). The fc*'* family contains local states \Si). gk 
indicates the number of different energies corresponding 
to this family. The k*^ family will be therefore repre- 
sented by gk states |m). We thus have J2k=i^k — ^4 
and X]fc=i5fc = 17- (For details, see TableP). 
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TABLE I: For each family, labels of different local states are 
given in column 2. Column 3 displays the generic form of 
local state : it stands for Nk states corresponding to a =] or 
j, and permutation of orbitals. gk and |m) are defined in the 
text. 



Because of the spin-flip term, off-diagonal elements of 
H\oc occur for local states belonging to families A: = 4, 7 
and 10. fc = 4 and /c = 10 families can be treated in 
a similar way. They can be described by six (2 x 2) 
blocks. As for the doubly degenerate case, singlet and 
triplet states have to be considered here. For instance, in 
the two-dimensional subspace formed by |>S'i4) = | ti 0) 
and IS'15) = I i,t,0), H\oc is represented by the (2 x 2) 
non-diagonal matrix 



U -J 
-J u 



Diagonalization leads to the first eigenstate IS'14) = 
(1514) + |<S'i5))/V2 which corresponds to £'14 = U — J 
and is the 5^ = component of a triplet (the corre- 
sponding — ±1 components are the states jcr, cr, 0) in 
the k = 3 family). The second eigenstate is the singlet 
state l^is) = (|5i4) - |S'i5))/V2 for which E15 = U + J. 
The notation is : |4) = IS'14) and |5) = \Si5). 

k = 7 family is composed by two (3 x 3) blocks of the 
following form : 

3C/-J -J -J 
-J 3C/-J -J 
-J -J 3U-J 



1526) = I T, T) and |527) = li, T, T), we have the eigen- 
states 



1^25) 
1^26) 
1^27) 
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V3' ''^ ' V3' 
with the corresponding energies 
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= 3U 
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The state |527) plays an important role for the dynami- 
cal properties of the triply degenerate Hubbard problem. 
Indeed this low energy state is one component of the max- 
imum spin multiplet 5 = 3/2. This fourfold multiplct is 
strongly stabilized by the Hund's exchange J. 

In table Hll energy eigenvalues E„i of |m) states are 
displayed. 
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TABLE II: Energy Er. 
on-site occupation. 



of each local state jm). Um is the total 



Considering the first block where \S; 



25/ 



T,T,i), 



H\oc is diagonal in the basis A, the effective Hamilto- 
nian Heft can therefore be solved by using the extended 
version of the non-crossing approximation presented in a 
previous work.* This approach has been used with suc- 
cess for various problems i^iSiLS 

Within NCA, propagators and self-energies for the 64 
local states \Si) are introduced. These propagators (and 
self-energies) are identical for similar local states. For 
instance | t,TjO)j I TjO, t) and | i,i,0) are described by 
the same propagator P3{uj) because of orbital and spin 
degeneracy. We have then only 17 different propagators 
Pm{i^)- However, all the 64 local states are needed to 



4 



express the Coaa and the c+j^ operators by Hubbard op- 
erators Xi j = \Si){Sj\. More precisely, we can write 

Coaa = ^D1■J\S^){SJ\ 

where is a (64 x 64) matrix and P"''^ is a (17 x 17) 
matrix. is calculated by developping the annihila- 

tion operator onto Hubbard operators and can be 
deduced by using the previously discussed degeneracy for 
local propagators Pm{uj). 

In the next section we present our results for the triply 
orbital degenerate Hubbard model. The density of states 
is obtained by solving the effective DMFT model by ap- 
plication of the NCA. 

III. RESULTS AND DISCUSSION 

In this section, we present the calculated densities of 
states for the triply orbital degenerate Hubbard model 
with a Hund's exchange invariant under spin rotation. 
The DMFT results are obtained for t = 1/V2 eV, U 
from 2 to 6 eV and for finite temperature T — 1000 K. 
We win take J « U/IQ or C//100. 

In order to understand the specific role played by 
Hund's coupling, and the nature of the different bands 
composing the density of states, we firstly investigate in 
the following sub-section the J = situation. The non 
vanishing J system will be studied in section Fill Bl 

A. No Hund's coupling 

Because of the triple band degeneracy, the electronic 
excitation spectra present already a rich structure for 
J = 0. This is shown in figure ^ where various concen- 
trations of charge carriers are investigated from n = 0.9 
up to n = 5.1. In this triply degenerate system, six elec- 
trons can in principle coexist at each site. For some plots 
(labelled by a star in the figure), the Fermi level position 
is indicated by a small circle. These results are in agree- 
ment with previous works on the MIT for degenerate 
correlated systems. In the strongly correlated regime, we 
find an insulating state for each integer filling. For fillings 
close to integer values, the narrow quasiparticle structure 
at the Fermi level is characteristic for a strongly corre- 
lated metal. With increasing chemical potential, an im- 
portant transfer of spectral weight is observed. The rel- 
ative band spectral weights can be interpreted in terms 
of probabilities. For instance, for n — 0.9 the excitation 
spectrum is mainly composed of two bands correspond- 
ing to the lower and the upper Hubbard bands (LHB 
and UHB). The number of occupied electronic states of 
the LHB gives the probability of removing an electron 
from the system, like in a photoemission experiment. 



The spectral weight of the UHB gives the probability 
of adding an electron in an empty state (inverse photoe- 
mission). For n close to 1 there is only one possibility 
for removing an electron, but there are five possibilities 
of adding an electron onto the same site. This is indi- 
cated in figurenfor various densities : (1, 5), (2, 4). . . For 
J = 0, the additional energy involved in this process is 
U , whatever the spin of the electron and the occupied 
orbital are. 




energy [eV] 

FIG. 1: Densities of states for various electronic occupations 
from n = 0.9 to n = 5.1 for the triply degenerate Hubbard 
model without Hund's coupling and for (7 = 6 eV. For some 
fillings (labelled by a star), the Fermi level is indicated by a 
small circle. 

Densities of states for n — 0.9,1.9,2.9,3.9 and 4.9 
are displayed in figure |21 with an energy shift ensuring 
that all Fermi energies are at the same point. Here, 
the transfer of spectral weight from the LHB to the 
UHB is obvious. The inset shows the low energy be- 
havior. We do not consider exactly integer values for 
n to prevent numerical problems when the density of 
states at the Fermi level is very close to zero. Nev- 
ertheless, our approach is numerically well defined for 
filling arbitrary close to integer values, with increasing 
numerical efforts. For all fillings, a quasiparticle peak 
is clearly present, but the effective mass of these quasi- 
particles is strongly dependent on the number of charge 
carriers. The correlated system with n = 0.9 is much 
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closer to the insulating state than the n = 4.9 situa- 
tion. The n = 0.9 effective mass can be obtained from 
the self-energy via ni* /m = (1 — dKeI]{uj) / duj)\^^Q and 
is approximately seven times the n = 4.9 effective mass. 
Actually, because of electron- hole symmetry, n = 1.1 and 
n — 4.9 give symmetric densities of states and there- 
fore the same effective masses for quasiparticles. Con- 
sequently, starting from n = 1, a 10% electron or hole 
doping leads to very different systems contrasting with 
the one band Hubbard model. For J = 0, we obtained 
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FIG. 3: Densities of states for various non-integer electronic 
densities with Hund's exchange and for (7 = 6 eV. (a) n = 
0.9, (b) n = 1.9 and (c) n — 2.9. Insets show low energy 
excitations. 



FIG. 2: Densities of states for various electronic occupations 
without Hund's coupling. For each plot U = 6 eV and the 
Fermi level corresponds to the origin of the graph. The inset 
shows the low energy behavior. 



a Mott metal insulator transition for each integer filling 
n — 1,2,3,4,5. This is consistent with previous theo- 
retical works using slave-boson formalism(2i2Si22i2& vari- 
ational method^ and quantum Monte-Carlo dynamical 
mean field theoryiiStH Nevertheless, the interesting ques- 
tion of the order of the transition cannot be investigated 
by the present approach because of the NCA pathology 
at low temperature^iiiS Indeed, the coexistence region 
extending to finite doping off half-filing observed for the 
doubly degenerate Hubbard model has been observed 
only at low temperature>i2i2S 

B. Non vanishing Hund's coupling 

Keeping a spin-rotational invariant Hund's exchange 
means that we have to handle explicitly the last term in 
expression (^. Recently, a very accurate study showing 
the importance of this term has been presented for the 
two-band Hubbard modeli^ Indeed, a completely differ- 
ent behavior occurs if one replaces the rotationally in- 
variant Hund exchange by an Ising-like one. 

Here we will focus on the specific effects of the Hund's 
exchange coupling on the electron dynamics for the 
strongly correlated metal and for the integer filling sys- 
tem which is insulating for large U. We will see that n 
plays a crucial role. 



In figure |2| the calculated densities of states are dis- 
played for various values ofn: n = 0.9, n = 1.9 and 
n — 2.9. They correspond to strongly correlated metals 
close to the insulating state. The narrow quasi-particle 
peak at the Fermi level is a signature of the large elec- 
tronic effective mass close to the MIT. For each n and for 
U = 6 eV, the calculation has been done for three values 
of Hund's parameter J = 0, U/100 and C//10. The (a) 
part of the figure shows that, for n = 0.9, the effect of 
a finite J value is restricted to high energy structures of 
the spectrum. This can be understood by considering the 
energies of the local states listed in tabled For n = 0.9, 
local states playing an important role in the electronic 
dynamics are those with occupations 0, 1 and 2. Larger 
occupations of local states are highly unlikely. The only 
J-dependent local energies are those for the two-fermion 
states 1 3), |4), |5) and |6). Transitions from singly occu- 
pied state to these states involve high energy excitations 
of the order of U. This J influence restricted to high en- 
ergy bands is specific to the singly occupied multi-band 
Hubbard model. Indeed, the physical situation is dra- 
matically different for n — 1.9 and n — 2.9 (part (b) and 
part (c) of figure 13). Here, a non zero J has a strong 
infiuence on high energy excitations as well as close to 
the Fermi level. For low energy excitations an impor- 
tant reduction of the spectral weight is observed. This is 
clearly visible in the figure insets and is in agreement with 
very recent results at half filling for the doubly degener- 
ate Hubbard model with Hund's exchangei^ Here, we 
found that this effect is also observed away from half fill- 
ing. Besides, this reduction is accompanied by the emer- 
gence of a multi-peak structure around the Fermi level 
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corresponding to the lifting of degeneracy of multiparti- 
cle states in table im This multi-peak effect is stronger 
for n = 2.9 which is consistent with the energy interval 
3 J between the three- fermion local states of spin 3/2 and 
1/2. If we approach n — the peak at around 3 J di- 
minishes such that we are left with a gap of the order of 
[7 at n = 3. 
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FIG. 4: Densities of states for various integer electronic occu- 
pations with Hund's exchange and for U = 2 eV. (a) n = 1, 
(b) n = 2 and (c) n — 3. 



The integer filling situation is examined in figure 01 
The stability of the metallic state is ensured by the in- 
termediate value [7 = 2 eV which is smaller than the 
critical Uc of the MIT for different integer fillings. J 
takes values and U/10. Here again, the effect of a non 
zero J is very weak for n = 1 (part (a)) and restricted to 
high energy excitations. A completely different behavior 
occurs for n — 2 and n = 3. For n = 2 (part (b)), the 
introduction of a non vanishing J = J7/10 is sufficient to 
destroy the quasiparticle coherent peak, which is consis- 
tent with the NRG results for the two-orbital Hubbard 
modeli^ This effect becomes stronger for the n = 3 sit- 
uation where a pseudo-gap is observed. Such a structure 
with a pseudo-gap at the Fermi energy is already present 
at half filling in the solution for the local impurity model 
proposed in ref.0. Here we show that this critical role of 
J, leading to a pseudo-gap appearance, can be observed 
away from half filling but the critical Coulomb and Hund 
parameters are strongly dependent on n. This property 
has important consequences on the nature of the differ- 
ent MIT in degenerate systems. A multi-band Hubbard 
model displays a Mott transition at all integer fillings. 
Nevertheless, by taking into account the Hund's exchange 
term, the different MIT present different low energy dy- 
namical behavior. It has been shown in previous studies 
that the critical correlation strengths Uc depend on n. 



For instance, for a multi-orbital Hubbard model with- 
out Hund's coupHng and for a ffat density of states, S. 
Florens and A. GeorgesS^ found that Uc is proportional 
to n(l — n). The maximum critical coupling Uc{n) for 
three orbitals is obtained for n = 3, the half-filled sit- 
uation, where orbital fluctuations are largest. This re- 
sult is consistent with previous work, in particular M.J. 
Rozenberg^° has found that Uc{2) > C/c(l) for the doubly 
degenerate orbital (with J = 0) Hubbard model by using 
Quantum Monte-Carlo calculations. 

Concerning the influence of the Hund's coupling J on 
the Mott transition at T = 0, an interesting generaliza- 
tion of the linearized dynamical mean field theory has 
been proposed recentlyi2& In addition to a reduction of 
Uc due to J, a qualitative change in the nature of the 
T = Mott transition has been observed for any finite 
J. The reduction of Uc has been given in details only 
for n = M where AI is the degeneracy. Here we find 
important differences between the half-filled case n = M 
and the n ^ M situation. Besides, our approach is able 
to explore dynamical properties by calculating the full 
densities of states. 




0.2 0.3 
J[eV] 

FIG. 5: Density of states at the Fermi level for different values 
of U from 2.4 to 3.6, as a function of J for n = 2. Note that 
the system becomes insulating for a critical J if [7 is large 
enough. 



Figure |31 shows the J influence on the density of states 
at the Fermi level for a filling n — 2 and for several 
U values. A non-zero J reduces low energy density of 
states. Interestingly, for moderate [/, there is a specific 
value of J, noted J™ for which the effect is the strongest. 
It is clear from figure [S] that increases with U and we 
have approximately Jm ~ U /lO. For large enough ?7, the 
increase of J leads to a vanishing spectral weight at the 
Fermi level. The resulting insulating phase is driven by 
Hund's exchange. Note that for n = 2, we found such an 
insulating phase for ~ 0.2 eV and only for U > 3.2 eV. 
A comparison with figure IHl shows the importance of the 
band filling n. 
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FIG. 6: Densities of states at £f as a function of J for various 
integer occupations n = 1, 2 or 3, and for U = 2.4 eV. Note 
that J plays an important role at half filling n = 3 but also 
away from half filling for n = 2. 



Indeed, in figure the dependence of the density of 
states at on J is plotted for n = 1, n = 2 and n — 3, 
and for a constant U = 2.4 eV. For n = 1, the influence of 
Hund's exchange is very weak. The only effect on the low 
energy excitations is a light increase of spectral weight. 
An appreciable decrease of the density of states at ep is 
found for n = 2 and U — 2.4 eV but the effect is much 
stronger for the half-filled system n — 3. Similarly to the 
n = 2 situation, the density of states at goes to zero 
for a critical value Jc ~ 0.2 eV but for n = 3 this takes 
place at a lower repulsion value U ~ 2.4 eV. 



IV. CONCLUSION 

We have proposed an approach based on the dynam- 
ical mean field theory which is able to describe impor- 



tant features of dynamical electronic properties of the 
triply degenerate Hubbard model. The self consistent lo- 
cal impurity problem is solved with a generalization of 
the non-crossing approximation. This impurity problem 
consists in a triply degenerate local correlated site (64 lo- 
cal states) embedded in an uncorrelated conduction band 
which is orbital dependent. As expected, for J = 0, we 
obtained a Mott metal insulator transition for each inte- 
ger filling namely n = 1,2,3, 4, 5. 

The influence of a Hund's exchange invariant under 
spin rotation has been studied in detail. In agreement 
with previous works, we found that a non-zero J implies 
a reduction of the low energy spectral weight. This effect 
is observed for integer as well as for non-integer fillings. 
In addition, we established the occurrence of multi-peak 
effects in the spectral function at the metallic side of the 
MIT. These effects give rise to structures with an energy 
scale J. An important result of our approach is the very 
strong dependence of all these effects on the filling. The 
destruction of low energy states is much more pronounced 
at half filling than at other fillings. Besides, for n = 1, 
this effect is not visible. Another interesting point is 
the occurrence of a particular value for which the 
reduction is maximum, J,„ increasing with U. 

It could be interesting to include in the present model 
different bandwidths for different orbitals to investigate 
the role of Hund's exchange in the so-called orbital selec- 
tive Mott transition^ Another possible perspective for 
this work is to understand how disorder in real mate- 
rials can affect the presented results concerning Hund's 
exchange influence. Finally, the pair hopping term be- 
tween different orbitals (orbital-flip term) can be studied 
in the framework of our approach and will be the subject 
of a forthcoming publication, as well as a more specific 
treatment of realistic spectra. 
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